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We study the arrest of three-dimensional flow in wet granular matter subject to a sinusoidal 
external force and a gravitational fleld conflning the flow in the vertical direction. The minimal 
strength of the external force that is required to keep the system in motion is determined by con- 
sidering the balance of injected and dissipated power. This provides a prediction whose excellent 
quality is demonstrated by a data collapse for an extensive set of event-driven molecular dynamics 
simulations where we varied the system size, particle number, the energy dissipated upon ruptur- 
ing capillary bridges, and the bridge length where rupture occurs. The three parameters of the 
theoretical prediction all lie within narrow margins of theoretical estimates. 

PACS numbers: 45.70.-n,68.08.Bc,05.70.Ln 



I. INTRODUCTION 

Sudden arrest of granular flows is an eminent prob- 
lem in the engineering sciences [l| as well as a challenge 
to the theoretical description of granular flows in a hy- 
drodynamic setting [2t9|. From the latter perspective 
it involves two challenges: (a) Appropriately incorporat- 
ing the role of dissipation arising from the particle in- 
teractions into the framework of the balance equations 
underlying hydrodynamic transport equations. And, (b) 
Addressing the role of shear stresses, of the spatial dis- 
tribution of stress, and of yield stress in systems where 
the flow is spatially anisotropic. 

For granular g^tems with purely repulsive interactions 
recent studies |10l - [l3t put severe constraints on hydro- 
dynamic descriptions of dense flows by pointing out a 
lack of scale separation of microscopic and relevant hy- 
drodynamic time and length scales. Among others this 
gives rise to a severe dependence of the effective mate- 
rial properties on the preparation history [l^ . In con- 
trast, hydrodynamic and continuum- mechanics consider- 
ations appear to provide a good description for granular 
systems where the hard-core collisions with restitution 
are augmented by (rev ersible) short-range attraction be- 
tween particles |15l-[l7l| . Arguably this is due to the sepa- 
ration of connectivity and rigidity percolation in response 
to attractive interactions This idealization of the 

particle interactions poj applies as long as high-impact 
particle collisions with high capillary numbers dominate 
the dynamics (see [2l|, [13 for recen t appli cations) . On 
the other hand, recent experimental j23l |24{ and numeri- 
cal |24l | work on slowly moving shear flow in dense gran- 
ular systems clearly underline the important impact of 
dissipation due to the hysteretic formation and break- 
ing of capillary bridges. Rather than accounting for the 
finite restitution in collisions and assuming reversible at- 
tractive forces, the present work therefore takes a com- 
plementary point of view: We explore slow flows in wet 
systems where dissipation is arising solely from the hys- 
teretic nature of the capillary interaction between the 
wetting liquid and the particles, i.e., it is due to the for- 



mation and rupturing of capillary bridges between parti- 
cles [m, [1^. The hard-core collisions are elastic. 

Shear forces that drive the flow can be modeled in vari- 
ous forms. Experimentally studying shear forces in gran- 
ular systems can be done, for instance, by constructing 
two counter-rotating cylindrical walls (see [1^ and ref- 
erences therein), by exploring a flow down an inclined 
plane [131 . Numerical models also addressed the shear 
flow induced by applying a cosine force field [1^ [28l - [30j . 
Similarly to the method of images, this may be used to 
implement zero flow velocity at the positions envisioned 
for the walls. 

Here, we will focus on the arrest of flow when the force 
driving the flow falls below a threshold value i^ex- In [lOl 
it was suggested that this critical force results from the 
power balance between the energy-injection rate result- 
ing from particles motion in the external force field, and 
the dissipation rate accounting for the rupture of cap- 
illary bridges between the particles. We will show that 
this approach might universally model the arrest in vastly 
different settings. To this end we describe the arrest of 
flow of wet, hard spheres in three dimensions where the 
motion in the third dimension is constrained by a hard 
wall at the bottom and a gravitational field in the ver- 
tical direction. For external driving forces close to i^ex 
all particles accumulate at the bottom of the cell such 
that the packing density lies only slightly above random 
close packing. In the following ways this setting is fun- 
damentally different from two-dimensional fiows at fixed 
density that were considered in Ref. [30j : 

• In the previously studied setting [s^l the particle 
density is prescribed by hard walls. Consequently, 
flow requires long-range rearrangements when the 
density is flxed close to random close packing. In 
that case the passing of particles in the flow re- 
quires a cooperative large scale rearrangement of 
large portions of the systems which is accompa- 
nied by the breaking of a large number of capillary 
bridges. 

• In contrast, the present study focuses on gravity- 
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FIG. 1. (a) Sketch of a system of width Lx = L and length Ly — 2L with periodic boundary conditions in the two lateral 
directions and an elastically reflecting top and bottom plate. The space-dependent cosine-shaped force field is indicated by 
the vertical band and by arrows. Additionally, a gravitational force, g, is acting downwards, (b) A snapshot of a simulation of 
2.15 X 10'' monodisperse spheres of diameter d that interact via a hysteretic square- well potential (see Sec. Ill A} . The cosine- 
shaped shear force field, Fs{x), is sketched on top of the figure. The color of each particle indicates its individual kinetic energy 
according to the color bar at the bottom left. The system size is L = 100 d and H — 12.5 d, and the filling height is h = 8.8 d 
for a shear force of amplitude -Fsmax = iOEch/L, i.e., -Fsmax = 0.15. 



confined granular beds of particles, where particles 
can pass each other via a slight, local expansion of 
the bed in the vertical direction. Only few capil- 
lary bridges must be ruptured in the resulting fixed- 
pressure setting, even though the density of the bed 
is always very close to random close packing. 

In spite of these considerable differences an informed im- 
plementation of the theoretical considerations developed 
in (soj provides an excellent description of arrest in both 
settings. 

We stress that the arrest of flow which is addressed 
here is conce ptu ally different from the jamming tran- 
sition [isl |3TI - I33j . Jamming is prevented here because 
the system can expand in the vertical direction where 
it is only bounded by gravity (see, for example 3J| and 
references therein). Furthermore, the transition is also 
distinct from those observed in earlier studies on cohe- 
sive granular materials, because we focus on the role 
of dissipation due to the hysteretic nature of capillary 
bridge ruptures, while previous studies [Tsl - flTj addressing 
the transition between fluidized and arrested states dealt 
with conservative attractive forces. They implemented 
dissipation via grain friction and a restitution coefficient 
smaller than one [l^ , or indirectly by treating the 
suspending fluid of attractive colloidal particles as an in- 
ert background [T5| . 

Our paper is organized as follows. In Sec. |TT] we de- 
scribe the system, our numerical approach, and the di- 
mensionless units adopted for modeling. Subsequently, in 
Sec.|lll]we work out the power injected into the systems 
by the external field, and the dissipated power due to the 
breaking of capillary bridges. The power balance pro- 
vides a prediction for the critical force i^ox- When larger 
forces are applied, energy injection dominates and there 



is sustained fiow [2^. Otherwise, the systems relaxes into 
an arrested state where the capillary bridges form a static 
network. In Sec. llVl the predicted parameter dependence 
of Fex is compared to the numerical findings. Based on 
only three scalar constants we can fully describe the de- 
pendence of the stability border in the four-dimensional 
parameter space spanned by the system size L, the parti- 
cle number N, the energy Ech dissipated upon rupturing 
capillary bridges, and the critical bridge length Scb where 
rupture occurs. (Here and in the following the subscript 
cb refers to capillary bridge.) The values of the three 
constants are interpreted and derived from the model in 
Sec.|Vl Finally, in the concluding Sec. IVII we summarize 
and interpret our main results. 



II. THE MODEL 

The system is confined in a rectangular cuboid of size 
L X 2L X H with periodic boundary conditions in x and 
y direction, and solid, refiecting walls in the z direction. 
This box contains N particles whose motion is confined 
in the z direction by a gravitational field of uniform ac- 
celeration, g = —gz, where z is the unit vector along 
the z-axis (see Fig.[T]). For the shear fiows studied in the 
present work particles never touch the upper wall of the 
container due to the gravitational confinement. 



A. Particle interactions 

In the present study we consider monodisperse spheres 
of diameter d, in order to suppress any additional dynam- 
ics arising from different particle sizes. After all, polydis- 
perse beads, which arc subject to small shearing forces. 
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segregate according to their size [28| . 

The dynamics in the simulation is calculated using a 
standard event-driven molecular-dynamics method which 
has been described in detail in Refs. [1^ [35l - l39| . For 
the sake of a self-contained exposition we only briefly 
summarize the particle interactions. Following [25| the 
effect of particle adhesion due to the capillary bridges is 
modeled as follows 

i. Capillary interaction gives rise to pair forces be- 
tween particles only. 

ii. When particles are not connected by a capillary 
bridge they feel no force when they approach each 
other. 

iii. They collide elastically, and upon collision a cap- 
illary bridge is formed instantaneously. When the 
particle separate, this bridge gives rise to an at- 
tractive force which is modeled as a potential with 
a depth Ech, and a finite range Scb- 

iv. The liquid bridge ruptures at the critical rupture 
separation, Scb, is modeled by removing the poten- 
tial well. When the particles approach the next 
time they feel no force again. In this manner the 
capillary bridge energy, Ech, is dissipated whenever 
a capillary bridge is removed (i.e., ruptures) after 
a collision. 

Work focusing on individual collisions [1^ liol - lisj reports 
a multitude of features for particle collisions involving 
capillary interactions that can not fully be captured by 
this model. On the other hand, the positions of phase 
boundaries of wet granular fluids appear to be universal 
in the sense that they only depend on Ech and Scb, and 
not on other details of the particle interaction [s^. For 
computational convenience all simulations shown in the 
present paper therefore adopt a hysteretic square- well po- 
tential, i.e., we use an event-driven algorithm where the 
potential takes the form of a square-well with hysteresis 
as outlined in i.-iv. 



B. Dimensionless units 

Masses are measured in units of the particle mass, to, 
distances in units of the particle diameter, d, and the time 
unit is fixed by measuring forces in terms of mg. Non- 
dimensionalised quantities are denoted by a hat. Unless 
stated otherwise the system size is L x 2L x = 60 x 
120 X 7.5, and the number of monodisperse particles is 
4.39 X 10*, resulting in a filling height ofh ~ 5.0. Further, 
the capillary interaction amounts to Ech = 3/8 and Sch = 
1/16. 
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FIG. 2. Time evolution of the kinetic energy, Ey, of motion 
parallel to Fs (upper set of lines), and the one in transverse 
direction, (lower set of lines) for a system of size L = 50 
driven by shear forces with amplitude Fsmax — 0.235 and 
0.245, respectively. 

C. Shear flow and arrest 

A shear flow is induced by applying a space-dependent 
external force field 

F,{x) = F,{x)y 

with Fs{x) = Fsniax ^OS (1) 

to the system [2^, HO, I46l - l48j , which accelerates particles 
along the y direction. The particles are initially homo- 
geneously distributed within the system with a Gaus- 
sian velocity distribution of mean granular temperature 
Tg/Ech = 40.0. 

For external forces with an amplitude, -Fsmax, slightly 
larger than the critical forcing, Fex, the system settles 
down into a stationary fluid flow whose local center-of- 
mass velocity follows the external field. In Fig. [IJb) we 
show a system with a relatively large filling height, h ~ 
8.8, and an external force only 11% above Fcx- In that 
situation shear bands form such that the capillary bridges 
in the region around Lx/4: and 3La;/4 are ruptured, while 
in the other parts the network of capillary bridges evolves 
only slowly. For shear forces smaller than Fcx the system 
eventually arrests in a solid state with a frozen network 
of nearest neighbors. 

D. Measuring Fcx 

The kinetic energy, Ey, of motion parallel to the driv- 
ing shear force will serve as the order parameter to dis- 
tinguish the dynamics [49| . In Fig.[2]its time evolution is 
shown together with the one of the kinetic energy, E^ , of 
the motion transverse to the external field. When started 



4 



in a state with high kinetic energy, the fluid first cools 
down in a manner closely reminiscent to free cooling [37| . 
Starting at i ~ 10 it has cooled so far that the accelera- 
tion due to the external field becomes noticeable over the 
initial kinetic energy. Depending on whether the exter- 
nal field is stronger or weaker than a sharp critical value, 
Fcx — 0.240, the fluid either gains sufficient energy to re- 
main in the fluid state forever (cf. (29j for the dynamics in 
that state), or it settles into the arrested state. The phase 
boundary, i^ex, is calculated as the mean value between 
the neighboring values of shear forces, Fismax, which ap- 
proach different states. For the data shown in Fig. [2] it 
amounts to i^ex — 0.240. 

The phase boundary hence corresponds to the smallest 
external force that still leads to sustained shear flow. We 
demonstrate in the following that the flow at this thresh- 
old corresponds to motion in a liquid state where the en- 
ergy injected by the external field is exactly balanced by 
dissipation due to rupture of capillary bridges when par- 
ticles move past each other. For all investigated systems 
the flow remains translationally invariant in y-direction; 
in accordance with the symmetry of the forcing. More- 
over, density of the liquid can not be distinguished from 
the random close packing density, 0rcp, and the height 
of the layer is spatially uniform. Consequently, the num- 
ber of particles, n{x), in a thin rectangular cuboid of size 
2L X h X dx aligned parallel to the external field takes 
a constant value, N/L, where N denotes the number of 
particles in the system and L the system size transverse 
to the flow. This finding is in line with the expectation 
that a fiow breaking the symmetries of the system or 
noticeably expanding the bed would give rise to higher 
dissipation. 



III. FLUIDISATION POINT OF SHEARED WET 
GRANULAR MATTER 



In this section we calculate the power, i.e., the energy 
injected per unit time into the kinetic energy of the par- 
ticles due to their motion in the external force field. In 
a steady state this power is balanced by the energy dis- 
sipation rate due to the inelastic particle collisions. For 
external forces close to i^ox the granular temperature is 
small. Therefore, effects due to the thermal motion may 
be neglected, and the energy dissipation rate is domi- 
nated by the rupture of capillary bridges when particles 
pass each other in the fiow field. 



A. Energy injection rate 

Let Vy{x) denote the flow profile of the flow es- 
tablished when applying the critical force Fs{x) = 

FsmaK cos(27ra;/L) = F^^ cos(27ra;/L). The input power 
that is injected by means of this external force acting on 



the particles is given by 

(-Pforcing) = / dx Vy{x) Fs{x) n{x) 



NF^ 
L 



/ dx Vy(x) cos —— . (2) 
Jo ^ 



Up to the factor N F^^ this amounts to the leading order 
Fourier coefficient of the expansion of the velocity profile 
Vy{x). Consequently, the injected power takes the from 

(^forcing) ^Cf N AVy Fe, (3) 

where Avy is the amplitude of the velocity field. Admis- 
sible values of C/ lie in the narrow range 0.5 < Cf < n/A. 
After all, the constant Cf takes the value 0.5 when the 
velocity profile is faithfully approximated by its first even 
harmonic, i.e., Vy(x) ~ Avy cos(27r x/L), and Cf = 7r/4 
in the other extreme case of plug fiow. 



B. Energy dissipation rate 

For each particle the creeping flow enforces a change 
of neighbors (in the direction of the fiow) with a rate 
7 = dvy/dx. Such a displacement goes along with an 
energy dissipation vEch due to rupturing v capillary 
bridges. The total power dissipated upon rupturing cap- 
illary bridges is thus given by 



(Pi 



bridge/ 



dx n{x) 



dv,, 



dx 



vE, 



cb 



(4) 



For every L-periodic function Vy{x) with a single maxi- 
mum this integral yields 



(Pbridgc) 



4 N Avy 
L 



cb 5 



(5) 



for spatially uniform, n{x) = N/L. 

For the present system rearrangements are achieved by 
a slight vertical expansion of the particle bed. Assuming 
that there is no height preference for the rearrangements, 
and that there are on average h/2d particles in the col- 
umn on top of the pair under consideration, a potential 
energy of 

U = Cu ^mg Ah 

is associated to the expansion. Here, Cjj accounts for 
the number of columns to be lifted, and Ah to the ex- 
pansion in height. Due to the frequent collisions in the 
dense bed the potential energy U is immediately dissi- 
pated into thermal degrees of freedom of the granular 
fiuid, and therefore it is not just a one-off investment. 
Rather, work has to be done against gravity each time 
particles move past each other. Multiplying this energy 
with the frequency of particle passages, ANAvy/L [as 
given by Eq. ([5|)], therefore provides a second contribu- 
tion to the energy dissipation, 

ANAvy ANAvy ^ Ah 
(Pgrav) = ^ ' U = ^ " Cu mgh . (6) 
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(-^forcing) — 


(-^bridge) ^" (-^grav) 




FcxL = 


a Ech + b mgh 


<^ 


FcxL = 


a Ech + bh . 
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with a = 








C'f 




b = 


2CuAh 









C. Predicting the critical force Fcx 

We now observe that in a steady state the total dis- 
sipation rate due to rupturing bridges and relaxing U 
into thermal degrees of freedom has to balance the input 
power. We hence obtain 



(7a) 
(7b) 

(7c) 



We hence predict that Fox is inversely proportional to 
the system size L, and that F^xL is a linear function of 
Ech and the filling height h. In the following section we 
compare these predictions to the results of the molecular- 
dynamics simulations. 



IV. PHASE DIAGRAMS 

In Fig. [3lja-d) we explore the dependence of Fox on 
the systems size, L, the filling height, h, the dissipated 
energy. Fob, and the rupture length, Scb- Plotting Fox 
as a function of the respective parameters provides sec- 
tions through the phase diagram: there is sustained flow 
for values of Eg larger than Fox, and flow is arrested for 
smaller external forces. 



A. Parameter dependence of Fcx 

(a) The variation of the system size L was done whilst 
keeping the aspect ratio of the container constant at 
Ly/Lj. = 2. At the same time the average particle 
number density was kept constant at (/) = 0.43 which 
means that the number of particles is changing in order 
to provide a fixed filling height, /i ~ 5. As expected 
from Eq. ([7]) the value Fex^ is constant. We find that 
FexF ^ 12 for fixed Fob = 0.375 and /i ~ 5. 

(b) The filling height, h, was varied by changing the 
number of particles in the system whilst keeping the ge- 
ometry of the simulation volume and particle interactions 
fixed. The filling height is estimated as the filling height 
in the solid state when assuming random close packing. 
As predicted by Eq. (O the dependence of Fox^ on h for 
a fixed Fob is linear. In the simulations we find 



the system size, L = 60, were kept constant. Fig. [3jc) 
shows the expected linear dependence of Fox-^ on Fob, 



FoxL ~ 6.0 + 1.2 h 



(8) 



FoxF ^ 24F, 



cb 



2.8. 



(9) 



for Fob = 0.375. 

(c) Whilst varying the capillary bridge energy. Fob , the 
rupture separation, Sob = 1/16, filling height, ft, ~ 5, and 



(d) The variation of the critical rupture separation, Sob, 
was done for a fixed aspect ratio of width vs. depth of the 
potential well, i.e. Fob was varied together with Sob at a 
fixed ratio of R = Fob/scb = 6.0. System size and filling 
height were fixed to F = 60 and h = 5, respectively. The 
resulting linear dependence 

Foxi ^ 16i?Sob + 6.0 (10) 

is shown in Fig. [3Kd). 

B. Consistency checks 

Merely finding linear dependences of Fox on h, Fob, 
and Sob is not sufficient to show that Eq. ([7]) faithfully 
describes the arrest of flow. One also has to verify that 
the linear functions arc mutually consistent for all sec- 
tions through the phase diagram. 

(e) Observing that the y-intercept in Eq. ^ may be 
written as 6.0 ~ 16 Fob for Fob = 0.375, and that the 
one in Eq. ([9]) is 2.8 ~ 0.56 ft for ft = 5, one realizes, 
that the linear dependencies, Eqs. ([5]) and can not 
be compatible with Eq. ([7]) at the same time. After all, 
0.56 disagrees with 1.2, and 24 disagrees with 16. In 
Fig. ^c) this is demonstrated by plotting F^xL versus 
Fob/ft: the data shown in Fig. [3fb) and Fig. ^c) lie on 
different straight lines. 

(f) On the other hand the data shown in panels (d) 
and (b) arc compatible. Eq. pop is commensurate with 
Eq. ([8]) since the latter has a y-intercept of 6.0 = 16Fob = 
16i?Sob with i? = 6.0 and Scb = 1/16, and since the y- 
intercept of Eq. ^ is 6.0 = 1.2ft with ft = 5.0. 

These findings suggest that Fob and ft arc not indepen- 
dent variables — as put forward by Eq. ([7]) when assum- 
ing that a and b take constant values. In the following 
we show that a consistent description of the numerical 
data can be achieved, however, by assuming that a has a 
weak, linear dependence on ft. It reflects that the num- 
ber of capillary bridges, v, ruptured in an exchange event 
depends on the filling height ft. 



V. DATA COLLAPSE 

Equation (fTU]) can only be consistent with Eq. © if 
its y-intcrcept comprises a contribution proportional to 
R, and if the slopes are adjusted accordingly, i.e., by 
decomposing the y-intcrcept as 6.0 = 1.2 ft ~ (0.55 + 
i?/10) ft with ft = 5 and i? = 6.0. Observing that R = 
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FIG. 3. The parameter dependence of the critical driving forces, -Fsmax, separating regions that lead to either solid or fluidized 
behavior. The symbols indicate the phase boundary between the two phases, and straight lines are the theoretical predictions 
described in the main text. For most of the data points the numerical error in the phase boundary is smaller than the size of the 
symbol. Panels (a)~(d) show the phase boundary between the solid and the fluidized state upon varying different parameters: 
(a) system size, L, (b) the filling height h, (c) dissipated energy, i?cb, and (d) critical rupture separation, Sch- In panel (e) 
Fe^L/h is plotted as a function of E^a/h in order to demonstrated that the data of panels (b) and (c) are not compatible with 
the prediction, Eq. ((Tjl. Panel (f) shows a data collapse of the data of panels (b) and (d), when Fe^L/h is plotted as a function 
Sch/h in order to demonstrate that Scb and h, rather than Ech and h may considered independent parameters to predict Fc^L. 
The method for performing the simulations is discussed in the text. 



-E'cb/scb this provides an improved prediction 



Determine 5 



Fc^L 



1 + X — a Ech 



with 5, ~ 16 , 
6-0.55 
X - 1/160 



(11a) 

(lib) 
(11c) 
(lid) 



which takes into account the increase of the number of 
bridge ruptures upon increasing h/scb due to the possibil- 
ity to rupture bridges in the vertically displaced column 
of particles on top of the site where a swapping event 
occurs. As demonstrated in Fig. |4] this prediction is in 
excellent agreement with all data. As a final step of data 
analysis we interpret the values of the fitting parameters 
a, b and x entering Eq. (|lip . 



In view of Eq. (j7b[) a is related to the number, v, 
capillary bridge ruptures in every particle exchange 



of 



4i/ 



(12) 



where we used the estimate Cf c± 0.5, as argued upon 
introducing this constant in Eq. ([Sj. Moreover, in order 
to take into account the correction term for the height 
dependence introduced in Eq. ([TI]) another factor 1 + 
/i/160scb must be accounted for. For the standard choice 
of parameters h — 5 and s^b = 1/16 we hence find 




16 



(13) 



Finally, in a granular bed that is expanding to allow par- 
ticles to pass each other the number v of rupture events 
must be larger than I, but still small. For = 3 we hence 
recover the value a = 16 reported in Eq. (jllbp . 
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(a) vary L 

(b) vary h 

(c) vary E^^^ 

(d) vary s^^ 




icb (/i-i + (160scb)-') 



FIG. 4. Master plot incorporating all data shown in the dif- 
ferent panels of Fig. |3] The dependence of the critical force 
i^ex on the system size, L, the filling height, ft, the energy dis- 
sipated upon rupturing capillary bridges, -Ecb, and the criti- 
cal rupture separation, Scb, is faithfully described by Eq. 
The solid line shows a theoretical prediction, Eq. involv- 
ing only three free parameters o, b, and x whose values can 
also explicitly be calculated. 



B. Determine b 

In order to see that the value for b ~ 0.55 is reason- 
able, too, we observe that subsequent rows of spheres in a 
elose-paeked structure are separated by a height distance 



'2/3 ~ 0.816 and that the saddle that has to be 
passed to roll from one minimum to a nearby minimum 
is at height /isaddic — 'n/S/S — 0.866. The minimum ex- 
pansion in height to move over the potential landscape 
set up by the layers below is therefore of the order of 
Aft, ~ /isaddie " ^cp — 0.05, and potential wells in a disor- 
dered, only slightly expanded random packing, will still 
be of the same order of magnitude. Based on this es- 
timate, on C/ ~ 0.5, and on Eq. ((7c)) the number of 
columns, Cjj , hfted a swapping event amounts to 



Cu ^ 0.55 



C 



f 



2Aft 



(14) 



Lifting a small number of columns lends plausibility to 
6-0.55. 



C. Determine x 

In order to gain insight into the order of magnitude 
of X we note that the lifted column needs to rupture 
bonds all along its walls. Hence, v is not merely depen- 
dent on its cross-section, as implied by assuming v to be 
independent of ft. The number of bonds that are bro- 
ken is then expected to scale linearly with column height 
and inversely proportional to Scb- After all, bonds are 



allowed to stretch to a finite length Scb, and the larger 
Scb the lower the likelihood that the dilation requires a 
given bond (along the vertical walls of the pile of parti- 
cles displaced vertically) to be broken. This leads to an 
additional number xvh/sch of bonds broken per column. 
As we saw above, it increases the number of bond rup- 
tures by about 50% which seems reasonable in a system 
with a packing height of only five layers and a rupture 
separation, Scb 1/16 — 0.065 that is larger than the 
height displacement. Aft = 0.05, of the column. After 
all, in such a situation only pre-stretched bonds are likely 
to rupture. Keeping this in mind, the small value of x 
may be understood as a result of incorporating a factor 
of Aft and a probability of about 1/4 to encounter a pre- 
stretched bond in a layer of a column that is vertically 
displaced. In this interpretation the average number of 
layers in the column amounts to ft/2. 



VI. SUMMARY AND OUTLOOK 

The present study substantiates the finding that the 
transition from a fluidized to an arrested state in wet 
granular matter arises when the dissipation rate due to 
the rupture of capillary bridges in the shear fiow can no 
longer be balanced by the power injection from the ex- 
ternal field. Earlier work [30| showed that this approach 
provides a comprehensive understanding of the transi- 
tion in the setting of a two-dimensional flow of bidisperse 
disks where the density is fixed by a confining box. The 
present work addressed the flow of a bed of monodis- 
perse spheres that are confined in the vertical direction 
by a gravitational field. Closely above the transition the 
fiow is reminiscent to a slow plastic flow of the bed in 
the direction of the applied field. The granular bed ex- 
pands only minimally in vertical direction. For the small 
fiow velocities gravity still keeps density, filling height, 
and pressure to values observed in an arrested packing. 
Also for this setting, which fundamentally differs from 
systems where the particles are confined by walls, the 
power balance provides an excellent prediction, Eq. (|lip , 
of the critical force where flow ceases. This is demon- 
strated by a master plot. Fig. 21 showing an excellent 
data collapse of numerical data obtained by varying four 
different characteristics of the system: the system size, 
the filling height, as well as the strength and critical rup- 
ture separation of capillary bridges. This data collapse 
suggests 

i. When the granular flow is confined by gravity, 
typically only 2-3 capillary bridges are ruptured 
upon swapping neighboring particles moving with 
slightly different speed in the direction of the ex- 
ternal forcing. This is a striking difference to 
wall-bounded flows [30j where this number diverges 
when the density approaches close packing. 

ii. In a gravity-confined granular bed the effortless 
passing of the particles is facilitated by a slight ex- 
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pansion of the granular bed where 2-3 columns of 
particles are lifted by a small amount to let the 
particles pass between neighboring potential wells 
arising from the corrugations formed by the layer 
below. The associated potential energy is also dis- 
sipated. 

iii. With a small probability additional capillary 
bridges are broken due to the slight expansions in 
vertical direction. 

The prediction of the flow-threshold involves only three 
constants characterizing the processes i.~iii. The values 
of these constants have been determined to a very good 
accuracy in Sec. |Vl 

The most remarkable finding of our study is that the 
critical force i^ox does not depend on the specific form of 
the flow profile. It can be calculated without specifying 
the hydrodynamic equations of the flow and determining 



their solution. We therefore conclude that the approach 
of balancing the energy input rate (due to the external 
force inducing the flow) and the dissipation rate (due the 
rupture of capillary bridges, when particles move past 
each other) provides a powerful framework to study the 
arrest of flow in wet granular materials where dissipation 
is dominated by capillary bridge ruptures. This approach 
provides a universal framework to predict the threshold 
for the arrest of flow, and it can be applied without need 
to determine flow profiles. Forthcoming work will explore 
this intriguing possibility also for flows in other geome- 
tries and due to other forcing. 
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